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Abstract 

A method is proposed for deriving dynamical equations for systems with both rigid and flexible com¬ 
ponents. During the derivation, each flexible component of the system is represented by a “surrogate 
element” which captures the response characteristics of that component and is easy to mathematically 
manipulate. The derivation proceeds essentially as if each surrogate element were a rigid body. Applica¬ 
tion of an extended form of Lagrange’s equation yields a set of simultaneous differential equations which 
can then be transformed to be the exact, partial differential equations for the original flexible system. This 
method’s use facilitates equation generation either by an analyst or through application of software-based 
symbolic manipulation. 
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Figure 1: Forces and moments generated at the bound¬ 
aries of a Bernoulli-Euler beam. 


1 Introduction 

Derivation of the exact partial differential equations de¬ 
scribing motion of systems with flexible elements is at 
best tedious and at worst so difficult as to drive the 
analyst toward other pursuits. A method is presented 
here which simplifies the process of equation develop¬ 
ment for a wide variety of system configurations. The 
key to the method involves modeling each flexible com¬ 
ponent in the system as a single dynamic element which 
captures the essential characteristics of the flexible mem¬ 
ber. This “surrogate element”, which can be treated 
substantially as if it were a rigid body, stands in for the 
flexible member through most of the derivation. Any of 
the standard continuum models may be used in the form 
of a surrogate element (i.e. Bernoulli-Euler beam, Tim¬ 
oshenko beam, etc.) The method also works for systems 
with rigid body rotations. [Singer 92] demonstrates the 
use of the surrogate element using Kane’s Method. 

First, the dynamic system is modeled with surro¬ 
gate elements substituted for flexible members. State 
variables are then assigned. Next, any of the conven¬ 
tional methods for deriving system equations can be ap¬ 
plied. When Lagrange’s method is employed, use of sur¬ 
rogate elements simplifies mathematical manipulations 
and eliminates the need to apply integration by parts. 
(Integration by parts is, generally, the most time con¬ 
suming step in the solution process for continuous sys¬ 
tems.) The equation generation process becomes very 
similar to that for systems made up of only rigid ele¬ 
ments. This representation and derivation strategy was 
developed to facilitate computer generation of equations 
of motion for flexible systems. Its use also provides ad¬ 
vantages when the equations are generated by the ana¬ 
lyst. 

2 The Surrogate Element 

A flexible component in a dynamical system will typ¬ 
ically exhibit the characteristics of both a mass and a 
spring. It will contain kinetic and potential energy and 
will exert forces and moments on adjacent components 
at its boundaries. If a flexible component is beamlike, it 
is frequently represented as a series of lumped elements, 
each with specified degrees of freedom. Bernoulli-Euler 
Beam elements and Timoshenko Beam elements are two 
common element representations. Each element of the 
series exerts a force and a moment only on its immediate 
neighbors. These forces and moments result from local 
deformations within the bodies. More specifically, the 
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Figure 2: This is an example of a surrogate element. 


load applied by an element depends only on the shape of 
the element at the point of application. (Of course, from 
a slightly different perspective, the shape is dependent 
only on the load.) An arbitrary number of similar ele¬ 
ments can be chosen to represent a flexible component. 

Just as the load applied by an element onto its neigh¬ 
boring element depends only on its shape at the bound¬ 
ary, the forces and moments applied by the entire flexible 
component onto adjacent lumped or flexible components 
are functions only of the flexible component’s shape (spa¬ 
tial derivatives) at its boundaries. For a Bernoulli-Euler 
Beam of length l, pictured in Figure 1, these forces and 
moments are [Crandall 82] 
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The sign convention chosen for these equations is as set 
by the directions of the vectors in Figure 1. 

The surrogate element is a model or representation of 
a flexible component which can be used instead of a tra¬ 
ditional continuous representation during equation gen¬ 
eration. The surrogate element exhibits the essential dy¬ 
namic characteristics of the component. It consists of a 
single elemental particle of the continuum (called a “rep¬ 
resentative particle”) and a spine defining the element’s 
shape. The element can exert forces and moments at 
the boundaries that connect it to the rest of the system. 
The surrogate element contributes a kinetic energy term 
because of motion of its elemental mass, and a potential 
energy term due to energy storage at the location of the 
representative particle. Because every internal particle 
of the continuum is topologically identical, the response 
of one particle can be used to represent that of the con¬ 
tinuum. The forces and moments at the boundaries en¬ 
force the connection of the surrogate element to the rest 
of the system. (Appendix A justifies this description of 
the surrogate element.) 

Generalized coordinates must be assigned to the sur¬ 
rogate element at each of the two boundary points and 



at the representative particle (Table 1 shows coordinate 
selection for a variety of surrogate elements.) Frequently 
the generalized coordinates which locate the ends of the 
surrogate element will also locate the bodies to which 
the corresponding beam is rigidly attached. 

To satisfy constraints set by Lagrange’s method, the 
generalized coordinates locating the representative par¬ 
ticle must be referenced to a coordinate system which 
is geometrically compatible with that of the bodies to 
which the surrogate element is connected. For example, 
in Figure 2: 
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The condition of geometric compatibility is discussed in 
Appendix B. Note that the surrogate element in Figure 2 
represents a Bernoulli-Euler beam and so the position of 
the representative point can be set with a single gen¬ 
eralized coordinate [Crandall 82]. If this were a surro¬ 
gate element for a Timoshenko beam, two generalized 
coordinates would be required to set both position and 
orientation at that point. 

While the Bernoulli-Euler beam has only one coor¬ 
dinate for describing the motion of each beam element, 
two coordinates are required at each boundary to specify 
position and orientation of the end of the beam. This 
fact is consistent with the continuous derivation of the 
Bernoulli-Euler beam model [Crandall 82] and is not a 
product of the surrogate element definition. (The surro¬ 
gate element concept simply highlights this anomaly.) 

For the case of the Bernoulli-Euler Beam, the surro¬ 
gate element has kinetic coenergy, potential energy, and 
Lagrangian given by 
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with terms as defined at the bottom of Table 1. 

Note that the kinetic energy term contains the ab¬ 
solute velocity of the surrogate element’s representative 
particle, v abs . If the neutral axis of the beam is sta¬ 
tionery, then v a b s = . However, if the beam is rotating 

or translating, the absolute velocity of the representa¬ 
tive particle must be determined (see Example 1). This 
is consistent with the specification of the kinetic energy 
of a lumped mass. Table 1 includes expressions for ki¬ 
netic and potential energy for several types of surrogate 
elements. 

In the derivation strategy developed below, surrogate 
elements are assumed to communicate with adjacent 
lumped components by applying loads to them. Table 1 
includes expressions for these loads as they should ap¬ 
pear in the equations. The lumped components will, in 


turn, communicate with the surrogate elements through 
the enforcement of boundary position and orientation 
compatibility. The loads applied by adjacent compo¬ 
nents on the surrogate element are captured implicitly 
by these geometric constraints. Consequently, they are 
not included explicitly in the derivation (see Appendix 
A). 


3 An Extended Form of Lagrange’s 
Equations 


The common expression of Lagrange’s equations: 
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3 = 1,2,..., n (9) 


is valid for Lagrangians of the form, 
L(€i,- ■■ ,£«; tsS-- •. Here E, (for j = 1,2,... ,n) 

are the generalized forces [Crandall 82]. When a system 
model includes surrogate elements, the potential energy 
terms become functions of the spatial derivatives of the 
C; §| and Appendix A demonstrates that a special 
form of Lagrange’s equations, 
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can be derived for systems whose models include rigid 
and surrogate elements. For this equation, L s is the La¬ 
grangian formulated to include surrogate element terms 
and the s Ej (for j = 1,2, ...,n) are the generalized 
forces which include the loads applied by the surro¬ 
gate elements onto adjacent components. The first two 
terms of Equation 10 are those of Equation 9. The 
last two terms are those required when spatial deriva¬ 
tives appear in the Lagrangian. This result is shown 
by [Lanczos 70] and [Meirovitch 90] and is also derived 
for use in Lagrangian density functions by Low and 
Vidyasagar [Low 88]. 


4 Derivation Strategy 

The steps in a strategy for using surrogate elements to 
derive exact partial differential equations of motion for 
systems with flexible components will be described here. 
The description assumes the use of Lagrange’s equation 
and is followed by a sample derivation chosen to illustrate 
application of the method. 

1. Model the system as a collection of rigid and flexi¬ 
ble elements. 

2. Temporarily substitute for each flexible component 
a surrogate element. 

3. Assign the required generalized coordinates to the 
surrogate elements. If the system configuration 
precludes motion in any direction, a generalized 
coordinate in that direction need not be included. 



4. Assign any additional generalized coordinates re¬ 
quired to specify the states of the system’s rigid 
elements. 1 

5. List the kinetic and potential energies of all surro¬ 
gate elements as presented in Table 1. 

6 . List the kinetic and potential energies of the 
lumped elements. 

7. From these terms for kinetic and potential energy, 
form a Lagrangian. 

8 . List the generalized forces. The portions of the gen¬ 
eralized forces applied by the surrogate elements 
to adjacent components are given in Table 1. 

9. Apply the extended form of Lagrange’s equations 
to the Lagrangian for each state variable. 

10. For those resulting equations for which the variable 
specifying the width of the representative particle 
is a common factor, divide through by it. In the 
example below, this variable is Ax. Note that these 
equations will be those for generalized coordinates 
representing displacements of representative parti¬ 
cles from their neutral positions. 

11. Integrate along the spine of each surrogate element 
those terms of the equations which still contain Ax, 
the variable specifying the width of the representa¬ 
tive particle of the surrogate element. Performing 
this step removes the surrogate elements from the 
analysis and yields the exact, partial differential 
equations of motion and the boundary conditions 
for the system as modeled in step 1 . 

5 Example One 

A flexible beam is fixed to a moving mass, mb, as shown 
in Figure 3a. A mass, m a , is connected to ground with a 
spring, mb is constrained to move only vertically and is 
connected with a spring and a damper to m a . The mass, 
m c is connected to ground through a third spring. A 
distributed force, f(x), is applied to the beam. System 
parameters are as presented in the figure. 

Figure 3b shows the system with a surrogate element 
substituted for the flexible beam. Because the left end of 
the surrogate element spine cannot rotate, a generalized 
coordinate is assigned for translation (< 72 ) but not rota¬ 
tion. Coordinates (74 and (75 represent the displacement 
and rotation of the right end of the beam respectively. 
Coordinate <71 locates mass, m a - Note that ( 72 , < 73 , and (74 
are all referenced to the same coordinate system; they 
all have been defined relative to mass, m a . 

The independent coordinates are chosen as: 

& = « i = 1,2,...,5 (11) 



a) Model of system with both rigid body and continuous 
elements. 



b) Model of system with surrogate element and assigned 
coordinates. 


Figure 3: System configuration for Example 1. 

x Be certain that the coordinates for the parts of the rigid 
system that connect to each surrogate element are defined to 
be compatible with the coordinates for the corresponding con¬ 
tinuum. Appendix B elaborates on this requirement. 




The kinetic coenergy, potential energy, and La- 
grangian are expressed as: 
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Note that V is now a function of (or ^Jr)- 

The new generalized forces, S H, ( , as defined in Sec¬ 
tion 3 are determined for each degree-of-freedonr: 
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For this example the terms to be substituted into La¬ 
grange’s equations (10) for £1 = qi are 
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The terms for £2 = 92 are 
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6 Conclusions 

A method has been proposed for deriving the partial 
differential equations of motion for systems of lumped 
and flexible components. When used in conjunction with 
Lagrange’s method, the process of equation generation 
is significantly simplified. Particularly, integration by 
parts is no longer required. Both system equations and 
expressions of boundary conditions result directly from 
application of this method. The method contains the 
same simple steps used to derive the equations of motion 
for lumped parameter systems. 

7 Appendix A 



In this section it will be shown that the application of 
the surrogate element method yields correct equations 
of motion and boundary conditions when applied to a 
dynamic system of rigid and flexible elements. First, 
equations will be derived for a general system using the 
variational formulation. Then the surrogate element for¬ 
mulation will be employed, yielding identical equations. 

A Lagrangian for a system of rigid and flexible com¬ 
ponents can be expressed as, 

n 

A = ^ [ Aflexiblej T brjg],] (37) 

i=i 

where Lflexibiej contains all of the kinetic and potential 
energy terms for the jth of n continuous elements, and 
Arigid contains the rest of the terms, those not dependent 
on spatial orientation within an element. Aflexibiej is the 
sum of the energy terms for the elemental particles in 
the jth continuum, so the Lagrangian can be written 

n 

i = EE LfjAzj + 

-^rigid (38) 

3 =1 Zj 

where the Lfj is a Lagrangian density function 

(see [Low 88] and [Meirovitch 90]) and Zj is a variable 
expressing the position of each element along the jth 
continuum. 

The components of the Lagrangian for the system are 
expressed in functional form as 

Af j = Lfj (ri,ri,ej(z j ),e j (z j ), 

e 'j( z j)’ e j{zj), t) (39) 

A r igid = A r jgid ( Tf , tj , 6j ( Zj a ), 6j (Zj ^), 

ej{ Z 3a)-> e j{ Z ja)i e j{ Z jb)i 

^( Z 3a)^'j( Z 3b)^) ■ ( 40 ) 

where the r* are state variables defining the positions of 
the m rigid elements (for i = 1, 2 ,..., m); the ej(zj) are 
coordinates which define deflections within the flexible 
elements; Zj a and Zj b are the locations of the endpoints 
of the jth continuous element; and ej{zj a ) and &j{zj b ) 




Table 1: This is a table of example surrogate elements 
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are the coordinates at those endpoints. Differentiation 
in time is denoted with a dot ( 6j(zj) = dej(zj)/dt); spa¬ 
tial differentiation is denoted with a prime (e'j(zj) = 
dej(zj)/dzj). The time derivative of e'j(zj) is denoted 
as &j(zj). 

The Lfj functions contain the terms in F which are to 
be spatially integrated, while the L r i g id function includes 
the remaining terms. The F r j g id function contains con¬ 
tinuous coordinates specified at the boundaries of the 
continuum (the ej(zj ) and ej(zj b )). For a system of 
only rigid elements, a single variable would be selected 
to specify the motion of each rigid component. However, 
for mixed systems, coordinates assigned to continuous el¬ 
ements specify the motion of some rigid components of 
the system so additional variables are not needed. For 
example, in Figure 3, ( 73 ( 0 ) = 52 ■ For the variational 
approach, ( 73 ( 0 ) appears in F r ; g id. With the surrogate 
element approach, qo would be used instead. 


7.1 The Variational Approach 


The summation of the energy terms in Equation 38 be¬ 
comes an integral when the continuous case is consid¬ 
ered, 
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The Variational Indicator may be formed from the 
Lagrangian, F and the work, W, done by generalized 
forces: 
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The S r i g id i are the generalized forces in the Sri directions 
(associated with the rigid body coordinates). The Efj 
are the generalized forces in the 8ej(zj) directions (asso¬ 
ciated with the flexible coordinates). The Sbcijj Sb C 2 j, 
Sbc 3 j, and “bc 4 j are the generalized forces in the 8ej(zj b ), 
8 e j(zj a ), 8e'j(zj b ), and <5e'(zj a ) directions respectively. 


These generalized forces are derived using the Principle 
of Virtual Work, and therefore do not contain internal 
force terms. The subscript Sb c i... Sbc 4 refer to the S’s 
associated with the four boundary coordinates of each 
flexible element. 

The variation operator is applied to each term, 
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Next, the variation is integrated by parts. The terms 
in the variation of time derivatives are integrated first in 
time. The terms in the variation of spatial derivatives 
are integrated first in Zj, 



d ( dL r igid ^ — 

~Jt \ W) rigidi 


+ 



3 



dL i] 

dej( z j) 


dzj 


8ri 


7 



d ( dL {j 
dt\de j(zj)J Zj 


9 ( dLt > \ A 

dzj \de'Azj)) ^ 


+ 


+ 


d 2 ( dL {j 
dzsAdeUzj) 


dzj 


3 tjdzj 


i\«i) 

Se j( z j) 

<9Eigid 


dt \ de-j ( Zj b ) J dej ( Zj b ) 


dL 


rigid 


d 


dL fj 


de'j(zjb) dz J \ de "( z n) 


j\*3b) 


I i ——'bcl j 


Se j( z i b ) 


E 


dL, 


rigid 


dL 


rigid 


3 

8L {j 


dt \de j ( Z j a ) J dej(z ja ) 


dLu 


de 'j( Z ia ) dz J \ de "( Z Ja) 


+1 


'bc2 j 


6e A Z ia ) 


E 


dL 


rigid 


dL 


rigid 


9L fj 

i ~ 777-7 + ^bc3 j 


de 'j( z jb ) 


dt ydefjizjt)J de'j(z jb ) 
6e j( z ib) 


E 


dL 


rigid 


dL 


rigid 


dLu 


de "( Z ia) 


dt \ d e'j( Z ia) J de 'j( Z ia ) 

dt 


+ 


- l bc4 j 


(44) 


Note that the summation in i is for each of the rigid 
body equation terms; the summation over j is for the 
terms due to the flexible components. Now the equations 
of motion can be extracted from Equation 44. The first 
set of equations arises from setting the coefficient of each 
br.i to zero, 


dLfj ^ z dL r igid 



dt \ dr, 


i = 1, 2,..., m (45) 


The next set of equations is derived from setting the 
coefficient of each Sej(zj) to zero. The summation over j 
disappears because each term of the summation contains 
a different 6ej(zj) for any given j , and therefore, becomes 
a separate equation. For the integral in Zj to be zero for 
arbitrary variations, the integrand must be zero yielding 
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The last four sets of equations are boundary condi¬ 
tions at the edges of the continuum corresponding to 
e j( z j a ) and ej(zj b ). Note that the summation over j is 
dropped for the same reason as in the previous equation. 
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In total there are m rigid body equations, n flexible com¬ 
ponent equations, and 4 n boundary equations. 

7.2 The Surrogate Element Method 

The equation derivation presented above is now per¬ 
formed with the same notation but using the surrogate 
element method outlined in this paper. The Lagrangian 
for this case is 

n 

£ - = EE Lf jL^Zj H- brinui (51) 

1=1 Zj 

with terms as defined in Equations 39 and 40. The La¬ 
grangian is formed from the equivalent rigid-particle sur¬ 
rogate element model. Equation 10 is applied to the 
Lagrangian, 
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The first m of the £* correspond to coordinates r, above; 
the next n of the £* correspond to coordinates e : j(zj); 
the last 4n of the £* correspond to coordinates ej(zj a ), 
e j( z j b ), e 'j( z j a ), and e'j(zj b ) respectively. Therefore, p = 
to + 5n, and j = k — m so that the Zj refer to the spatial 
variables along the continuum corresponding to £*. The 
are the generalized forces defined in terms of the 
generalized forces from the variational approach given 

by 


-rigid* 
-'f (k—m) 


-'bcl (k—m—n 
H - Stable/c 


—hc2 (k—m—2n) 
^ +' = 'tablefc 


- J bc3 ( k—m—3n) 
+ 1 Stable k 


-‘hcA(k—m—An) 
Stable k 


for k = 

1, 2, . .. , TO 
for k 
■ + 1 , 

to + 2 ,..., to + n 
for k = 
m + n+ 1 , 
to + n + 2 , 

..., to + 2n 
for A: = 

m + 2n+l, (53) 

to + 2?i + 2, 

..., to + 3n 
for k = 
to + 3?r + 1, 
to + 3?i + 2, 

..., to + 4n 
for k = 
to + 4?i + 1, 
to + 4?r + 2, 

..., to + 5n 


where the S ta bi e * are obtained from Table 1 and can only 
be non-zero for coordinates at the boundaries of the con- 
tinua. 

The £* correspond to three different variable types; 
the rigid body coordinates, the continuous coordinates, 
and the continuous coordinates at the boundaries. The 
evaluation of Equation 52 is treated separately for each 
case. For k = 1,2,..., to, Equation 52 becomes, 


dLtj . 0F, 

-rr^Azj -- 

orj or. 


EE 

j 

d (dLt 


rigid 


dt \ dri 


Azj 


j A*; 

d ( <5L ri gid 
dt \ dfi 


■"rigidj 


i = 1,2,..., to (54) 


The third and fourth terms in Equation 52 are zero be¬ 
cause the £* are rigid coordinates. 

For k = toT 1, to+2, ..., m+n , Equation 52 becomes, 


<9Ffj d ( dL {j \ 

dej(zj) ^ + dtydejizj)) 

d ( dLu \ A a 2 / OLu \ 

+ dz 3 [de'^zj)) AZj dzf yde'j(zj)) 

= EfjAzj j = 1,2,... ,n . (55) 


For k = to + n +1, to. + n + 2,..., to + 5n, Equation 52 


becomes, 


d f dL ri gi d \ 
dt \dej(zj b )) 

d { F F r jgi,i \ 

dt XdejizjJJ 

d j F Frigid j 

dt ydefajh)J 

d ( dL tigid 

Vde'feJ 


^Frigid _ 

^ e j( Z jb ) 

“bcl j + Fi j 

SFjigid 

5e i(*J«) 

“bc2j + F 2j 

5F rig i(j 

de 'Mj b ) 

“bc3y + FF| j 

^Frigid 

“bc4j + 


J = 1,2,... ,n (56) 

J = 1,2, • • • ,n (57) 

J = l,2,...,n (58) 

j = l,2,...,n (59,) 


where the F\j, F 2 j, M\j, and M 2j are obtained from 
Table 1. The Table defines these values as the forces and 
moments acting from the continuum onto the adjacent 
components. The values used in creating the table are 
given by 


Fi 


F 2 


Mi 


M 2 


dLfj d f dLfj \ 

de'j{z jb ) dzj \de'](z jb ) J 

dLfj d f dLfj \ 

de 'j( z ia) + 9z j \ de "( z j a )) 
dLfj 

de "( z Jb) 

dLfj 

■ 


(60) 

(61) 

(62) 

(63) 


Following the procedure presented in Section 4, A Zj is 
canceled from every term of Equation 55 to yield Equa¬ 
tion 46. Next, all terms in Azj that remain in any of the 
equations are integrated along the Zj. These terms are 
all found in Equation 54. The integration yields Equa¬ 
tion 45. Finally, Equations 56-59 can be rewritten with 
substitutions from Table 1 for the Fi, F 2 , Mi, and M 2 
to obtain Equations 47-50. 


8 Appendix B 

The surrogate element derivation above assumed that 
the coordinate describing the location of the representa¬ 
tive particle in a continuum be geometrically compatible 
with the coordinates for adjacent elements at the bound¬ 
aries of that continuum. Geometric compatibility means 
that the variation of the surrogate element coordinates 
and the variation of the boundary coordinates be equal. 
Consider the example shown in Figure 4. Coordinate r 2 
at x = 0 equals r\, and r 2 at x = l equals r%. However, 
q -2 at x = 0 is equal to zero, and q 3 equals the sum of <71 
plus q 2 at x = l. The method outlined above assumes 
that the consistent coordinate set (ri, r 2 , 7-3) be chosen. 
The requirement on the coordinate selection, 
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= dr 1 


(64) 




Figure 4: Example using coordinates with different refer¬ 
ence frames. The j*i, V 2 , r$ frame is referenced to ground. 
The gi and (73 are referenced to ground while (72 is refer¬ 
enced to the left hand mass. 

and 

i = Sr 3 (65) 

is satisfied. The force terms on the end of the contin¬ 
uum are multiplied by the variations of the coordinates 
at those boundaries (< 5 r 2| x=0 and Sr 2 \ x=l in this exam¬ 
ple). The force terms acting on the rigid components 
are multiplied by the variations of the rigid component 
coordinates (^rq and Sr 3 in this example). The method 
outlined above implicitly equates these two sets of forces. 
However, if the two sets of variations are not equal for 
the chosen coordinate system, then the equating of these 
force terms can not be justified. 

The requirement discussed above can be relaxed by 
realizing that the force terms can be corrected when us¬ 
ing a selection of coordinates that is not geometrically 
compatible. In the example above using the 91 , <721 93 
system, Sq 2 \ x=(j = 0. The force term F\, which would 
be applied to the mass on the left when using a geo¬ 
metrically compatible coordinate system, is not added 
to s Si. Also, Sq 2 \ x= i = 693 — Sqi. Therefore, for La¬ 
grange’s method to provide correct equations of motion, 
F 2 , the force at the right end of the beam must be added 
to the generalized force term for the <73 equation ( s Ss), 
and subtracted from the generalized force term for the 
<71 equation ( s Si). This modification is convenient when 
applying the surrogate element method to more complex 
systems such as those with rotating flexible components. 
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